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Abstract 



We describe a new mechanism for the generation of the baryon asymme- 
try of the universe during a first order electroweak phase transition. The 
mechanism requires the existence of two (or more) baryon number carrying 
scalar fields with masses and CP violating mixing which vary with the Higgs 
field expectation value. This mechanism can be implemented using squarks 
in supersymmetric theories or using leptoquarks. Our central observation is 
that reflection of these scalars from a bubble wall can yield a signiflcant net 
baryon number flux into the symmetric phase, balanced by a flux of oppo- 
site sign into the broken phase. For generic parameter choices, scalars with 
incident energies in a specific, but not narrow, range yield order one refiec- 
tion asymmetries (between the probability of refiection of the scalars and of 
their antiparticles). The interesting energies are those for which there are two 
propagating scalars in the symmetric phase but only one in the broken phase. 
Electroweak sphaleron processes drive the baryon number in the symmetric 
phase toward zero, but do not act in the broken phase. Our estimate of the 
resulting baryon asymmetry is consistent with cosmological observations for 
a range of mass parameters and CP violating phases in a supersymmetric 
implementation, as long as the bubble walls are not too fast and not too 
thick. 
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I. INTRODUCTION 



Understanding how the baryon asymmetry of the universe (BAU) could have arisen from 
matter-antimatter symmetric initial conditions remains a central challenge at the interface 
between high energy physics and cosmology. The observed light element abundances, com- 
bined with the theory of big bang nucleosynthesis, yield the constraint that the ratio of the 
cosmological baryon number density to the cosmological entropy density lies in the range 
3 X 10^^^ < ub/s < 9x 10^^^ |l|]. As has been understood since the work of Sakharov 
0, a theory that explains the generation of a net baryon asymmetry must include depar- 
tures from thermal equilibrium, and must include processes that violate baryon number 
conservation, charge conjugation symmetry C, and CP symmetry, the product of charge 
conjugation and parity symmetry P. These conditions can be satisfied by electroweak pro- 
cesses in the early universe, raising the possibility of an electroweak explanation of the 
BAU[3-28]. In particular, Sakharov's requirements are all satisfied during a first order elec- 
troweak phase transition [7-26]. As the universe cools through the transition temperature Tc, 
non-equilibrium conditions exist near the bubble walls that are sweeping through the uni- 
verse, converting the high temperature "symmetric" phase to the low temperature "broken" 
phase.Q Baryon number is violated in the electroweak theory as a consequence of the anomaly 
PP| in processes in which the sphaleron barrier is traversed . The rate for baryon number 
violating processes is exponentially small at zero temperature and in the broken phase at 
the critical temperature T^, but not in the symmetric phase [^,3T|. C is maximally violated 
in the electroweak theory, and CP violation enters via the Cabibbo-Kobayashi-Maskawa 
(CKM) matrix. 

Meeting Sakharov's conditions, as the minimal standard model does, is necessary for 
generating a cosmological baryon asymmetry, but need not be sufficient for generating an 



^By convention, throughout the rest of this paper we refer to the phase at T = Tc that is smoothly 
connected to the equihbrium phase at T > Tc as the symmetric phase, and the phase at T = Tc 
that is smoothly connected to the equihbrium phase at T < Tc as the broken phase, even though 
there is in fact no distinction between the symmetry of the two phases. 



asymmetry as large as that observed. Indeed, in order to explain the observed BAU, it 
appears that one needs to extend the standard model for at least two reasons. First, although 
mechanisms using CP violation from the CKM matrix alone have been constructed [1^, they 



seem to be unable to generate the observed BAU WM. Hence, CP violation beyond that in 
the standard model seems to be required. Second, regardless of the mechanism by which 
baryon number is generated, if it is to survive after the electroweak phase transition, baryon 
number violating processes must be sufficiently suppressed for T < Tf.. This imposes a lower 
bound on the expectation value of the Higgs field in the broken phase at T = |^ which 
is not satisfied in the standard model with an experimentally allowed Higgs mass. This 
bound can, however, be satisfied in various extensions of the standard model, including the 
minimal supersymmetric standard model with one top squark lighter than the top quark 
p3| ; throughout this paper, we will assume that the theory is suitably augmented such that 
the transition is strongly enough first order that the bound is satisfied. Since the energies 
that are potentially relevant to electroweak baryogenesis will be increasingly explored in the 
coming decade by present and planned accelerators, it is of great interest to explore new 
weak-scale physics that can explain the observed BAU. 

In this paper, we introduce a new mechanism for generating the BAU at the electroweak 
phase transition. Our mechanism requires augmenting the standard model by the addition 
of (at least) two baryon number carrying complex scalar fields (pi and 02 with masses of 
0{Tc). The two fields must be coupled by off-diagonal terms in their (Hermitian) mass- 
squared matrix which include a CP violating phase. In this way, we introduce CP 
violation beyond that in the CKM matrix. Furthermore, we require that depend upon 
the Higgs field expectation values in the theory so that the mass eigenvalues and eigenstates 
are different in the symmetric and broken phases. The requirements just sketched can be 
implemented in a variety of extensions of the standard model. Perhaps the most appealing 
possibility is that 0i and 02 are the SU (2) singlet and SU (2) doublet top squarks in a 
supersymmetric extension of the standard model. Another possibility, of interest in light 
of the recent HERA anomaly P2[, is that 0i and 02 may be weak-scale leptoquarks whose 



masses receive contributions from couplings to the Higgs field. Throughout this paper, 
we focus on the baryogenesis mechanism rather than on model building. For definiteness, 
however, we present our mechanism and results taking the 0's to be squarks,^ and defer 
discussion of other possibilities to the concluding section. 

In a supersymmetric theory, the mass-squared matrix depends upon fi(T) and V2{T), 
the temperature dependent expectation values of the two Higgs doublets Hi and H2 that 
give mass to the down and up type quarks, respectively. The off-diagonal terms in are in 
fact zero in the symmetric phase, where f 1 = f 2 = 0. In the broken phase, the off-diagonal 
terms are complex and therefore CP violating. More generally, we define 0i and 02 as the 
eigenvectors of in the symmetric phase, and note that the mass eigenstates in the wall 
and in the broken phase are linear combinations thereof. 

We will be interested in the reflection and transmission probabilities of 0's incident on the 
wall from the symmetric phase. Because of the CP violating phases in M^, the probability 
R12 for an incident 0i to be reflected back into the symmetric phase as a 02 is not the same 
as Rj2, the probability for an incident antiparticle 0i to be reflected as a 02. We will show 
that this reflection asymmetry AR = R12 — R-j2 results in a net flux of baryon number 
from the bubble wall into the symmetric phase, compensated by a flux of the opposite sign 
into the broken phase. Note that the measure of CP violation in the model is the spatial 
variation of the phase of the off-diagonal terms in M^. A spatially constant phase can be 
rotated away by a spacetime- independent unitary transformation on the 0^; hence the phase 
must vary spatially if AR is to be nonzero. 

The central observation of this paper is that the reflection asymmetry AR can be large, 
approaching 1, over a broad range of incident energies, if the phase of the off-diagonal term in 

changes by (9(1) as the bubble wall is traversed from the symmetric phase into the broken 
phase. As already noted, because of the dependence of on vi and V2, the eigenvalues of 



^Previous treatments of electroweak baryogenesis in supersymmetric theories include those of 



Refs. |ll,19,23-2e]. 



vary within the bubble wall. We assume that the larger of the two eigenvalues in the 
broken phase is greater than the larger of the two eigenvalues in the symmetric phase. This 
means that there is in general a range of incident energies E such that in the symmetric 
phase, both eigenvalues of are below E"^, while in the broken phase, there is only one 
eigenvalue below E'^. Therefore, there are two propagating modes with energy E in the 
symmetric phase and only one in the broken phase. Consider a 0i incident upon the wall 
from the symmetric phase with an energy in this range. As it begins to penetrate the 
wall, it evolves into a linear combination of the position-dependent eigenstates of the matrix 
M^. Since only one mode can propagate in the broken phase, there is a position within 
the wall at which one mode is totally reflected. The reflected mode, upon re-emerging into 
the symmetric phase, is some linear combination of (pi and 02 which includes a significant 
02 component if there is significant mixing. Another way of understanding what is special 
about the range of energy under discussion is that at each energy in this range, there is 
one linear combination of incident 0i and 02 that is totally reflected. For this reason, both 
Ri2 and Rj2 are generically large, and AR is also large unless the CP violating phase is 
small. We will refer to this range of incident energies as the "enhanced reflection zone."0 
By comparison, for incident energies above the enhanced reflection zone, for which there 
are two propagating modes in both the symmetric and broken phases and throughout the 
bubble wall, we find that AR is nonzero but is generically many orders of magnitude smaller 
than one. The width in energy of the enhanced reflection zone is comparable to the amount 
by which the masses change between the two phases; in the example we present, the width 
of the enhanced reflection zone is 90 GeV. In Section II, we present the parametrization of 

appearing in a supersymmetric theory. We then set up the calculation of AR, leaving a 
detailed presentation of the method of calculation to the Appendix. We evaluate AR and 
explore its dependence on parameters in M^. 



We will show that for energies in the enhanced reflection zone, scalars incident upon the wall 
from the the broken phase do not yield a reflection asymmetry and hence do not contribute to the 
BAU. 
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Standard model quarks also have an enhanced reflection zone, as discussed by Farrar and 
Shaposhnikov , although its width is only of order the strange quark mass. The resulting 



BAU is small [jT^, essentially because the light quarks have mean free paths much shorter 



than their Compton wavelengths. We defer to Section IV a discussion of the suppression 
due to the finite mean free path of the heavy scalars we employ in our mechanism; the 
suppression is not severe. 

In Section III, we integrate AR against the appropriate thermal distributions for incident 
0i's and 02's to obtain the baryon number flux injected into the symmetric phase. If the wall 
velocity is zero, or if the masses of 0i and 02 in the symmetric phase are equal, we find 
that the baryon number flux due to incident 0i's is cancelled by that due to incident </'2's. 
As long as Vw ^ and the masses in the symmetric phase are not degenerate, we obtain 
a nonzero baryon number flux. The larger the fraction of the thermal distributions for 
incident 0's lying in the enhanced reflection zone, the larger the baryon number flux will be. 
The final element in the mechanism involves electroweak baryon number violating processes. 
These drive the baryon number density in the symmetric phase toward zero. Because they 
do not act in the broken phase, the final result is a net baryon asymmetry of the universe 
whose magnitude we estimate in Section IV. Our mechanism yields a BAU consistent with 
observation if the scalars have nondegenerate masses of order Tc in the symmetric phase 
and if the bubble walls are sufficiently thin and slow. We discuss open questions and model 
implementations in Section V, and note there that an enhanced reflection zone can arise in 
leptoquark models, and thus is not peculiar to supersymmetric theories. 

To close this introduction, we contrast our mechanism with the charge transport mech- 
anism, pioneered by Cohen, Kaplan and Nelson |7;12], further developed by many authors. 



and used to estimate the BAU generated during the electroweak phase transition in super- 
symmetric theories[19, 23-26]. Our mechanism can be seen as a modification of the charge 
transport mechanism. We make explicit comparisons with the results of Huet and Nelson 
||T9| obtained using the charge transport mechanism and find that our mechanism can yield 
an n^/s consistent with observation for smaller CP violating phases. The central difference 



is that in our mechanism, we generate a flux of baryon number into the symmetric phase, 
whereas in the charge transport mechanism a flux of another quantum number, often left- 
handed baryon number minus right-handed baryon number, is generated. This axial baryon 
number can be washed out by QCD processes before it has time to bias electroweak baryon 
number violating processes Because our mechanism generates a baryon number flux, 
it is immune to QCD interference of this kind.0 This contrast is particularly germane in 
light of the recent demonstration that the rates for the relevant QCD processes are sig- 
nificantly larger than previously expected |3^. Various authors have noted the possibility 



that a baryon number flux may be generated, but this has always been assumed to be a 
small effect. This is in fact true for incident energies such that the number of propagating 
modes is the same on both sides of the bubble wall. For a generic mass-squared matrix M^, 
however, there is a broad region of incident energies in which fewer modes propagate in the 
broken phase. We observe that this leads to large reflection asymmetries, and consequently 
to a large baryon number flux into the symmetric phase, yielding an efficient mechanism for 
generating a BAU consistent with cosmological observations during the electroweak phase 
transition. We will refer to the electroweak baryogenesis mechanism we propose as the scalar 
baryon number transport mechanism. 

II. M^, AR, AND THE ENHANCED REFLECTION ZONE 

We begin this section by presenting the parametrization of appropriate when (pi and 
(j)2 are left- and right-handed top squarks and set the stage for the calculation of AR. We 
then describe the dependence of AR upon the incident energy and upon parameters in M^. 



^We should note that there are scenarios in which baryogenesis via the generation of an axial 
baryon number current can be immunized against suppression due to strong sphalerons. One 
example |jl^ requires that the left- and right-handed top squarks and the left-handed bottom 
squark have symmetric phase masses comparable to the temperature while the other squarks are 
heavier. Another example ||3^ requires the formation of a squark condensate just above Tc- 
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As discussed in the introduction, the CP violation that we exploit results from the mixing 
between two scalars as they traverse the bubble wall separating regions of symmetric and 
broken phase at T = T^. In this background, terms in the potential that couple the baryon 
number carrying scalars 0j to the Higgs fields Hi and H2 give the (pi spacetime-dependent 
masses and mixings which can be encoded in a 2 x 2 Hermitian mass matrix. In a super- 
symmetric theory, the scalars 0j are squarks. Top squarks (stops) are the most promising 
candidates to play the role we envision for the 0i's because they can be light without violating 
experimental upper bounds on neutron and electron electric dipole moments (EDMs). In- 
deed, Cohen, Kaplan, and Nelson have recently advocated supersymmetric models in which 
CP violating phases are 0{1), but observable EDMs do not arise because the first and sec- 
ond generation squarks have masses in the tens of TeV . Although it is conceivable that 
the scalar baryon number transport mechanism could be implemented using first or second 
generation squarks, it seems likely that stops will yield the largest contribution to the BAU. 
The stop mass matrix is fSSl 

/ m^^ + mf + m| cos 2(3 (^| sin^ Ow) mt{Ae''^^ + yue^'^^ cot 13) \ 

y mtiAe-'^f^ + ne-'^fB cot (3) mf^ + mf + m| cos 2/3 - | sin^ Ow) / 

Here, rrit = XtV2 is the top quark mass, m| = g'^{vl+V2)/2 is the Z-boson mass squared, rhtR 
is the soft SUSY-breaking mass for the SU (2) singlet stop, rhtL is the soft SUSY-breaking 
mass for the SU{2) doublet stop, tan/5 = V2/V1 is the ratio of the Higgs field expectation 
values, /ie*''''^ is the (complex) mass in the Higgs potential coupling the two Higgs fields, and 
Ae*'^-* is a complex soft SUSY-breaking term. We will see that the scalar baryon number 
transport mechanism works best if rhtL and rhtR are both 0{Tc) and differ by about 10-30%. 
Mass differences of this order can arise due to renormalization group evolution down from 



some high energy scale at which rritL = rritR. Indeed, in the models of Ref. ||36|, rutL and 
rhtR differ by 20%. 

During most of the existence of the expanding bubble, its wall can be treated as flat 
because its radius of curvature is much larger than its thickness, so the mass matrix will 



depend only upon one spatial direction, which we take to be the x direction. We choose 
the convention that the region of large negative x is the symmetric phase and that of large 
positive X is the broken phase. Both vi and V2 vary across the bubble wall. A complete 
calculation of these profiles is beyond the scope of this paper, although a treatment using 
the resummed one-loop temperature-dependent effective potential is possible [^. As is 
conventional, we make a simple choice in terms of a single width parameter, hoping that 



this captures the essential physics. Following Ref. |2J, we choose profiles such that the Vi{x) 
are x- independent for x < —w/2 and x > w /2 and are sinusoids for —w/2 < x < w /2. That 
is, we define the profile function 

x< -w/2 

p{x) = { I + i sin (^) -w/2 <x<w/2 (2.2) 



and then for vi, we take 



2 ' 2 

1 ' ' x>w/2 



vi{x)=v',p{x) , (2.3) 

which varies smoothly from zero in the symmetric phase at a; — oo to fj* in the broken 
phase at a; — > oo. The parameter w characterizes the width of the wall separating the 
two phases. Although in reality the profiles behave exponentially for |a;| ^ w, most of 
the interesting physics happens where x is varying most rapidly and p{x) is as good a 
parametrization as any. The choice ( |2.2| ) is convenient numerically, as it allows us to impose 
boundary conditions at x = ±w/2 rather than at larger 

It is crucial to our mechanism that V2/V1 vary across the bubble wall so that the overall 
phase of the off-diagonal term in is not constant. We take 

V2{x) = f i(a;) tan/5(x) 

= vi{x) [p{x) tan p{x)) tan(/3 - A(3)\ (2.4) 

so that tan/3(a;) varies from tan(/3 — A(5) in the symmetric phase to tan/5 in the broken 
phase. We will take tan/3 = 2 in our estimates. Our results do not depend sensitively 
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on this choice. The appropriate choice for w° in ( |2.3|) is not (250 GeV) cos^, the value it 
takes at T = 0, but rather the value it takes in the broken phase at T = T^. This can 
be calculated as a function of parameters in specific models, but we will simply use the 
reasonable estimate f° = (2/3)(250 GeV) cos/?. (Our results do not depend sensitively on 
the choice of prefactor.) In order for the baryon asymmetry generated (by any mechanism) 
during the electroweak phase transition not to be wiped out, (f?)^ + (^2)^ must be larger 



than . Estimates for A/5 exist in specific models and range from 0.01 — 0.03 p^ , |38[| to 
0.25 [^, but there are certainly no experimental constraints on this parameter. In order 
for the overall phase of the off-diagonal terms in to vary with x, that is, in order for 
to introduce CP violating effects, we must have 

A/5 7^0 

/^^ = ^A-^B^^ ■ (2.5) 



EDM experiments may constrain A(y9 in some models |3^, but it has recently been noted 



that in other models is constrained to be small while if a is essentially unconstrained. 
Any constraints on the CP violating phases are weakened if the first and second generation 
squarks are heavy. Finally, note that (pA can be generation-dependent. Hence, there is no 
model-independent constraint on /\ip for third generation squarks. 

As a concrete example which will serve as a visual aid for much of our subsequent 
discussion, in Figure 1, we plot the eigenvalues of as a function of x for the parameter 
set rhtL = 110 GeV, mtR = 90 GeV, A = /i = 100 GeV, tan/? = 2, tan(/3 - A/5) = 1, 
that is A/? = 0.32, Ay? = n/2, niz = 91 GeV, = 175 GeV, and sin^ 9w = 0.23. With 
these parameters, the zero temperature masses of the two top squarks are 141 GeV and 
243 GeV. In our example, we take the wall width to be w = (4 GeV)^^. Although the 
critical temperature and the wall velocity f ^ play no role in the calculation of AR, we 
mention for completeness that when they enter in the next section, we will use = 100 GeV 
and = 0.1 as representative values. It is conventional to write the wall width in terms 
of the temperature, that is, w = 25/Tc in our example. Henceforth, we write as T when 

9 



200 




60 I 

-w/2 xo +w/2 

Position 

FIG. 1. Mass eigenvalues as a function of position. The two curves show the square roots of the 
eigenvalues of M^, with parameters given in the text, as a function of x. w is the wall width. Hsi 
and iJ,s2 are the mass eigenvalues in the symmetric phase, and //^i and 111,2 are the mass eigenvalues 
in the broken phase. For any incident energy in the shaded range fis2 < E < 111,2 , there is a position 
Xo at which the upper eigenvalue crosses E. 

it causes no ambiguity. The final relevant parameter is the (p mean free path I, which will 
appear in Section IV where we estimate I ~ 10/T ~ 0.4^. This completes the enumeration 
of the parameters specifying the "canonical" example for which we will quote results in the 
following. We will, of course, describe the effects of varying each of these parameters at 
appropriate points in the discussion. 

We wish to follow a (p particle from the thermal ensemble in the symmetric phase that 
last scattered somewhere away from the wall and is now incident upon the wall. Implicit 
in this scenario is the assumption that the mean free path of the scalars in the plasma is 
long compared to the wall width w. This assumption is likely false, but we nevertheless 
use it in this section and the next, deferring our treatment of the suppression due to finite 
mean free paths for the particles to Section IV. The particle impinges upon the wall and 
is reflected or transmitted, and then resumes its thermal motion in the plasma on one side 
of the wall or the other. During the time between the last scattering before reflection or 
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transmission and the first after, tlie particle propagates freely, feeling only the changing 
expectation values of the Higgs fields which are encoded in the mass matrix. We calculate 
the reflection coefficients and their asymmetry in the rest frame of the wall; we will boost 
the resulting baryon number flux to the plasma frame when we calculate it in Section III. In 
the wall frame, energy is conserved upon traversing the wall since the mass matrix is time- 
independent. The reflection coefficients can be calculated by solving the time independent 
Klein-Gordon equations 

0,(x) = O (2.6) 

because the time dependence of solutions is simply an overall exp{iEt). In general, we will 
take a basis in 0i and 02 such that far from the wall in the symmetric phase, is diagonal. 
This has already been accomplished, since m^^ and m^^ are the only terms in ( |2.1| ) which are 
nonzero in the symmetric phase. When calculating the baryon number flux in Section III, 
we will consider particles incident upon the wall with momenta that are not perpendicular 
to the wall. However, their reflection coefficients will depend only upon the component of 
their momentum that is perpendicular to the wall; hence, it suffices to compute reflection 
coefficients for normal incidence. Therefore, the problem of finding reflection coefficients 
is effectively a one-dimensional scattering problem, and equations (|2.6|) become ordinary 
differential equations for (t>j{x). 

Consider the mass matrix whose eigenvalues are shown in Figure 1. The behavior of 
the reflection coefficients is qualitatively different for particles incident from the symmetric 
phase with energies E < ^s2, 1^32 < E < fib2, and E > fib2- For E < /is2, there is only one 
propagating mode in both the symmetric and the broken phases. In general, we denote the 
reflection coefficient for a (pi from the symmetric phase reflected back into the symmetric 
phase as a (pj by Rij for particles, and denote the corresponding quantity for antiparticles by 
Rjj. For E < /is2, however, the only reflection coefficients are -Rn and Rjj;. Since the CPT 
conjugate of the reflection of 0i to 0i is the reflection of 0i to 1, we see that Rn — ^tTi 
and there is no reflection asymmetry. Before proceeding to higher energies, note that for 
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a different mass matrix it may be tfie case tfiat > //s2- In tfiis case, for E < there 
are two propagating modes in the symmetric phase and none in the broken phase, so both 
01 and 4>2 must be totally reflected. We now argue that in this circumstance, the reflection 
coefficients again cannot be CP violating. Unitarity implies that for total reflection, 

Ru + i?i2 = 1 (2.7) 
^TT+^T2 = 1 (2.8) 

We see that unitarity, together with Ru — Rj^, implies that R12 — R12 and hence there is no 
reflection asymmetry. For the rest of this paper, it is implicit that references to /is2 should 
be replaced by references to /i^i if the mass matrix is such that /ibi > A*s2- We have shown 
that particles incident from the symmetric phase with energy E < yield no reflection 
asymmetry. 

Now consider incident energies iis2 < E < 11^2, for which there are two propagating 
modes in the symmetric phase and only one in the broken phase. As shown in Figure 1, for 
any energy in this shaded range there is a point Xq at which the larger eigenvalue of 
equals E"^. This means that there is a particular linear combination of incident 4>i and 4>2 
that evolves by mixing as it propagates through the wall in just such a way that upon arrival 
at Xq it is purely in the mass eigenstate with eigenvalue E^, and is therefore totally reflected. 
Since one linear combination of 0i and 02 is totally reflected by the wall, both R12 and R21 
are large in this energy range, given sufficient mixing. In order to obtain large asymmetries 
between the reflection coefficients for particles and for antiparticles, the individual reflection 
coefficients must of course be large, making this zone of enhanced reflection a promising 
place to look for large AR. Without CP violation, the same linear combination of incident 
01 and incident 02 mixes to become the mode with eigenvalue E'^ at Xq, and is totally 
reflected. This implies that Rj2 = R12 and 7?2T = R^i- However, if the off-diagonal term in 
has a spatially varying phase, then the linear combination of 0's that is totally reflected 
is different from that for the 0's, and we expect AR ^ 0, as we conflrm explicitly below. 

Before proceeding, note that R^i — R21 by CPT and therefore 
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Energy (GeV) 

FIG. 2. AR as a function of incident energy E, for two different wall widths w. The enhanced 
reflection zone is shaded. 

AR = Ri2 -Rv2 = Ri2 - R21 . (2.9) 

This simplifies calculations by giving the reflection asymmetry in terms of reflection coef- 
ficients for particles only. Another simplification is that in the enhanced reflection zone, 
AR is the only possible asymmetry because there is no asymmetry due to particles incident 
from the broken phase. Denoting the single propagating mode in the broken phase by 03, 
since R32 = R^ by CPT, unitarity requires that incident ^s's cannot yield baryon number 
asymmetric transmission into the symmetric phase. We give a careful explanation of the 
calculation of AR in the Appendix; henceforth in this section, we focus solely on the results. 

In Figure 2, we plot AR vs. incident energy E for two values of the wall width, w = 5/T 
and w — 25/T, with all other parameters as in our canonical example. The enhanced 
reflection zone is apparent, with AR ~ 0.1 for w — 5/T. Even greater values of AR 
are obtained for A/3 and A(p that are larger than our canonical values. (For example, 
for A/3 = 1 and A(p = 2.5, we flnd AR ~ 0.3.) We have therefore confirmed that in the 
enhanced reflection zone, large reflection coefficients yield large reflection asymmetries in the 
presence of C(l) CP violating phases. For incident energies below the enhanced reflection 
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FIG. 3. AR as a function of wall width w for incident energy E = 150 GeV, which is within 
the enhanced reflection zone. 

zone, AR = 0. At higher energies, that is for E > ^1^2, there are two propagating modes in 
both the symmetric and the broken phases and throughout the wall. In this regime, AR ^ 0, 
but it is extremely small. For w = 5/T, AR = 0.0047 at E ^ 200 GeV; AR = 2.5 x 10"^ 
ai E = 250 GeV; and AR = 1.3 x 10~^ at E = 300 GeV. For these energies, transmission 
coefficients are very close to 1 and reflection coefficients (and therefore AR) are very small. 

For energies above iib2, unlike for those in the enhanced reflection zone, there is the 
possibility of a nonzero asymmetry arising from particles and antiparticles incident from the 
broken phase and transmitted through the wall into the symmetric phase. This asymmetry 
is very small, as we now argue. Just as the Rij are very small in this energy range, the 
reflection coefficients for reflection of particles incident from the broken phase back into the 
broken phase are very small also. Since any asymmetry associated with transmission into 
the symmetric phase must be balanced by an asymmetry in reflection, we conclude that 
even though the transmission coefficients are ~ 1, their asymmetry is very small. 

We now discuss the dependence of AR on the parameters in the problem, beginning 
with the wall width w. Note that in order to obtain a nonzero AR, there must be some 
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region in x in which modes incident from the symmetric phase can mix, and hence feel the 
effects of the CP violating terms, before arriving at xq where one mode is totally reflected. 
This implies that AR = for an infinitesimally thin (step function) wall. In the thin wall 
limit, Ri2 and Rj2 are large in the enhanced reflection zone, but they are equal. In Figure 
3, we plot AR at E = 150 GeV versus the wall width w for our canonical M^. We see that 
AR ~ w for w small relative to the inverse mass scales in the problem. At large w, the 
number of wavelengths per wall width grows and therefore R^ and Rj2 decrease, and so 
does AR. The asymmetry AR peaks at about w = 5/T, one of the values we have chosen 
to plot in Figure 2. This is an unphysically thin wall — estimates for w range from 10/T 
to 100/T. (The physics determining w is presented, for example, in Refs. [^,^].) In our 



canonical example we follow Ref. pj] and use w = 25/T, and we have plotted AR for this 
wall width in Figure 2. If w is in fact smaller than 25/T, our final result for the BAU is 
enhanced, while for thicker walls, it is somewhat suppressed. 

Let us now consider the effects of varying the mass parameters in M^. We define m and 
Am through 

Am 

fJ'si = nitR = m — 

Am 

fis2 = mtL = m + . (2.10) 

For simplicity, we will always take fi = A and tan /5 = 2. (The optimal choice for (3 should 
be such that A n cot /?. We find, however, that for A = varying P from 0.5 to 4 changes 
our results by at most 10%, so the dependence on /? is not significant.) We have investigated 
the dependence of our results on m. Am and A. Of course, we do not vary the m| and 
m^ terms in M^, and these should be thought of as setting the energy scale. Varying m 
between 50 and 150 GeV while holding Am/m and A fixed changes AR by less than 20%. 
Increasing m further leads to a suppression of AR because it increases all the eigenvalues 
relative to 1/w. For A = 0, there is no mixing between 0i and 02, and AR = 0. Increasing 
A from to 200 GeV holding all else fixed yields a monotonically increasing AR, but in 
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going from A = 100 GeV to ^ = 200 GeV, the increase in AR is less than 10%. Of the three 
mass parameters we have varied, Am has the biggest effect. AR is maximized for Am — 0. 
However, we will see in the next section that the baryon number flux is identically zero for 
Am = 0. As Am is increased, AR falls; by a few percent at Am/m = 0.2; by about a factor 
of two at Am/m = 0.8. We defer a plot displaying the dependence of our results on Am to 
the next section. 

Finally, we come to the CP violating phase. We have verified that for small A(p and A/3, 
AR is linear in both quantities. As we have noted, there is no model-independent constraint 
on Aip; therefore, in Figure 4 we show AR over the entire range of A(/9. We see that AR is 
approximately linear in A(p for \A(p\ < 2. We find that AR is linear in A/3 over the range 
—0.2 < A/3 < 0.5. The values we have been using in our canonical example — Acp — 7r/2 
and A/3 = 0.32 — are within the linear regime. For convenience, in subsequent sections, 
we often quote results assuming that A/3 and Aip are in their linear regimes, as is likely the 
case for A/3, but not necessarily for A(p. Fortified by our understanding of how AR in the 
enhanced reflection zone varies with wall width, masses, and phases, we are ready for the 
next step in our computation of the BAU. 

III. FROM ASYMMETRIC REFLECTION COEFFICIENTS TO BARYON 

NUMBER FLUX 

We now compute the baryon number flux into the symmetric phase by integrating AR 
against the flux densities of the incident particles. The flux, which we denote by J^b, receives 
two contributions. One contribution, Fg^s, is generated by the asymmetry in the reflection 
of particles and antiparticles incident on the wall from the symmetric phase and is associated 
with AR. The other contribution, Fj,^s, is caused by the asymmetry in the transmission 
of the particles and their antiparticles incident on the wall from the broken phase. We 
have seen, however, that this asymmetry (unlike AR) is zero for energies in the enhanced 
reflection zone, and is small (like AR) at higher energies. We therefore neglect F^^s and 
obtain 
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FIG. 4. Ai? as a function of for incident energy at the center of the enhanced reflection 
zone. 

J'b ^ \f,^, . (3.1) 

The factor of 1/3 arises because the particles have baryon number 1/3. 

Before deriving an expression for Fg^s, we make some general observations. The pro- 
cesses we are studying are invariant under time translations (although not under time re- 
versal!). Therefore, energy is conserved upon reflection from or transmission through the 
wall. We denote the components of a particle's momentum in the rest frame of the wall that 
are perpendicular and parallel to the wall by pi± and pi\\, respectively. Recalling that is 
diagonal in the symmetric phase, we have 

= + p2 1 + = + p2 1 + , (3.2) 

where Pi\\ = \pi\\\- Since the processes of reflection from and transmission through the wall 
are invariant under spatial translations in the directions parallel to the wall, pi\\ is conserved. 
For example, for a (pi reflecting into a 02, Pi\\ = P2\\- Equation ( |3.2| ) then yields 

pI± + /^si = pI± + /^s2 = £± > (3-3) 
17 




defining a conserved quantity e±. We now see how to translate results obtained in the 
(1 + l)-dimensional treatment of the previous section into the full (3 + l)-dimensional setting 
appropriate here. What was called E in previous sections is in fact e±. AR depends on e±, 
and the enhanced reflection zone is given by ^s2 < e± < Hh2- 

Before the arrival of the wall, the 0's are in thermal equilibrium with momenta distributed 
according to the Bose-Einstein distribution. This equilibrium distribution defines a rest 
frame, the plasma frame, which is different from the rest frame of the wall. In order to 
compute J^B, we need the flux density of (pi particles (equivalently, (pi antiparticles) incident 
upon the wall with a given momentum in the wall frame. Throughout this section, we 
continue to assume that the mean free path of particles is larger than the wall width, 
deferring the discussion of the effects of the falseness of this assumption to Section IV. 
Hence, the incident flux density we require is given simply by 

?>PijE 



ME,Pi±) 



(3.4) 



where 7 = l/^/l — v"^. The factor of 3 appears because there are 0's with each of 3 different 
colors in thermal equilibrium. The argument of the exponential arises because, as mentioned 
before, the particles are initially in thermal equilibrium in the plasma frame, not in the wall 
frame. To get Fg^si we must integrate A_R = R12 — Rj^ = R12 — R21 against /i over the 
three-momentum of the incident (pi and integrate R21 — R2T = R21 — R12 = — Ai? against /2 
over the three-momentum of the incident 02 and add up the results. Thus, we obtain 
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(27r)2 7p2±min 27r 

where PiXmin is the Pix_ such that e_[_ = fis2- Note that Fg^g is the flux seen in the plasma 
frame. Transforming from the wall frame back to the plasma frame yields the overall factor 
of 1/7 in ( p.5|) . Upon performing the integration over pjii, we obtain 
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(3.6) 



where we have used (|3.3| ) in the form e±de± = pi±dpi± = P2±dp2±- 

From (|3^ ), we deduce that Fg-^s = for v^] = 0. This is to be expected, since baryo- 



genesis requires out of equihbrium conditions, and hence requires Vw ^ 0. (Note that one 
can derive an expression for Ff,^s similar to that for Fs_,s and show that F^^g = when 
Vu, = 0.) The vanishing of J^b with can be more directly understood as follows. We 
present the argument in 1 + 1 dimensions; the generalization to 3 + 1 is trivial. We have 
shown that R12 — Rj2 = —{R21 — -^21)- For 1;^ = 0, the number of (pi particles incident upon 
the wall with an energy E greater than /i<j2 is equal to the number of incident 02 particles 
with the same energy. (There are of course 0i particles with /x^i < E < /is2, but they do not 
yield a reflection asymmetry.) Therefore, for v^, = 0, the contribution to J-'b due to incident 
01 's is exactly cancelled by that due to incident 02 's. Now consider a moving wall, where 
Vw 0. 0j particles incident upon the wall with a given energy E in the wall frame have 
energy {E — VwPi±) in the plasma frame, in which they are initially in a thermal equilibrium 
distribution. If fi^i 7^ fis2 then pi± 7^ p2± for a given E, and the number of 0i and 02 
particles with incident energy E in the wall frame is not the same. We see that in order to 
upset the cancellation between the contribution to J-'b due to incident 0i's and that due to 
02 's, we need both f «, 7^ and figi 7^ P-s2- The asymmetry in the reflection coefficients can 
only yield an asymmetry in the baryon number flux if the wall is moving and if the scalars 
are not degenerate in mass in the symmetric phase. 

As discussed in the previous section, AR depends on parameters in and on the 
wall width. The flux J-'b depends on these parameters as well as on the wall velocity and 
temperature. We will be interested in the regime in which and Am/m are small compared 
to 1. This may seem surprising, given that we have just argued that jFg is zero for Vy^ ^ or 
Am/m 0. The reason is that, as we have seen in the previous section, AR is a decreasing 
function of Am and, as we will see in the next section, the BAU is proportional to J-'b/v'^ 
for realistic wall velocities. In order to gain intuition about (p.6|) it is useful to pretend that 
AR is energy independent for /is2 < e± < fib2, and then expand in v^, and Am/m to first 
order in both, obtaining 
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FIG. 5. Temperature dependence of the baryon number flux. In this plot, T varies and 
w = 25/T for ah T. We keep aU parameters except T and w fixed. We plot Tb/T^ vs. m/T 
to facilitate comparison with (|3.7| ). 



1 Am f m\'^ , 
v,„ ( — ) Ai? 
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ydy 



(3.7) 



T3 47r2 - m v r/ Jrr./T [ev -l]^y^ - {m/ry 

We perform all our calculations using (|3.6|) , not the expansion (|3.7|) , but the expansion is 
useful for understanding the qualitative dependence of J^b on the parameters. 
With all parameters as in our canonical example, ( p.6|) yields 

^ ~ 9 X 10-^ , (3.8) 

and with w = 5/T instead of 25/T, we obtain a result which is a factor of two larger. We 
have verified that J^b is linear in to within a few percent for < 0.6. In the regime in 
which jFg is linear in A/?, Aip, and v^, (p.8|) becomes 



:^|~2xlO-S^A/5A(^. 



(3.9) 



Turning now to the temperature dependence, in Figure 5, we plot Tb/T^ vs. m/T, 
varying T and keeping all parameters in the mass matrix fixed. Since Ai? does not depend 
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FIG. 6. Dependence of the baryon number flux on Am/m, with Am varying and all other 
parameters fixed. 

on T, we can partially understand this plot by noting that in ( pj.?] ) jFg/T^ ~ (m/T)^ at 
small m/T and ~ exp(— m/T) at large m/T. This does not completely describe the Figure, 
however, because as we vary T, we have kept w = 25/T; this means that w changes with 
respect to the parameters in the mass matrix. We conclude from Figure 5 that the BAU 
generated by the scalar baryon number transport mechanism is largest for m ~ T. 

In Figure 6, we plot jFg/T^ vs. Am/m. It is linear in Am/m for small Am/m and falls 
at large Am/m because, as we noted in the previous section, AR decreases with increasing 
Am/m. We see that using Am/m = 0.2 as in our canonical example yields a reasonable 
estimate for AR over the range 0.1 < Am/m < 0.3, but for Am/m outside this range, 
the BAU is suppressed. For the scalar baryon transport mechanism to be efficient, we need 
scalars with symmetric phase masses of 0{Tc) that differ by 10-30%. 

IV. ESTIMATING THE BARYON NUMBER OF THE UNIVERSE 

In the preceding sections, we have shown how to calculate the baryon number flux J-'b 
carried by particles that is injected into the symmetric phase by the motion of the bubble 
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wall. To this point, we have described the quantitative solution of a well-posed problem. 
Given a mass matrix, a critical temperature, a wall profile, a wall velocity, and making the 
assumption that the mean free path is long compared to the wall width, a quantitative 
calculation of Tb is attainable. In this section, we sketch a qualitative estimate of the 
cosmological baryon to entropy ratio, n^/s, that results from the flux jFg. Our treatment 
is admittedly crude and can be improved, for example along the lines of that of Huet and 
Nelson [|l^, but we leave this for future work. We organize the estimate of the final result 
as follows. First, we estimate the mean free path / and the suppression of jFg that results 
from the finiteness of l/w. Then, we estimate the scalar baryon number density (baryon 
number in the form of 0's) that results from jFg. This in turn leads to a quark baryon 
number density which biases the electroweak sphaleron processes acting in the symmetric 
phase, resulting in a net baryon asymmetry of the universe. 

Before the wall arrives, the (/)'s in the symmetric phase are diffusing about with a mean 
free path in the x direction which we call I and a mean velocity in the x direction between 
scatterings which we call v^. The one-dimensional diffusion constant is defined as 

D = lv^. (4.1) 

The mean one-dimensional velocity is ~ 0.7 for particles with m = T. Joyce, Prokopec, 
and Turok have estimated that in the symmetric phase, the diffusion constant for rela- 



tivistic strongly interacting particles is D ~ 6/T. Most of the scatterings contributing to D 
are interactions with gluons in the plasma, and D is inversely proportional to the relevant 
interaction cross-section which decreases with increasing m. This suggests that for m = T, 
the diffusion constant is somewhat greater than 6/T, but we leave a complete calculation 
for future work. In this paper, we simply use the conservative estimate 

Note that an increase in m/T leads to an increase in I, and potentially to a larger final 
result. Of course, in increasing m/T, one pays a penalty in the exponential suppression in 
(p. 71), and we will see below that making the 0's nonrelativistic exacts other costs as well. 
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FIG. 7. Dependence of the baryon number flux on l/w, with I varying and ah other parameters 
fixed. 

We now give a crude estimate of the suppression due to the finiteness of l/w. As noted 
in Section II, estimates for w range from 10/T to 100/T. In our canonical example, we use 
w = 25/T, yielding l/w ~ 0.4. The essence of the effect is that when (p particles reach xq, 
the point in the wall at which one mode is totally reflected, they have only been travelling 
(and mixing) freely for a distance of order /. To incorporate this, we redo the calculation 
of AR as follows. For each incident energy, we find xq and choose as incident states the 
mass eigenstates a distance / to the left of Xq. These then propagate only a distance / before 
reflecting, and so experience less CP violating mixing than in the case where I is infinite. 
The result is a suppression in AR at each energy. We in fact find that this suppression is 
rather energy dependent, being larger for energies close to fib2- In evaluating jFg, therefore, 
we must re-evaluate the integral ( p.6| ). The result is shown in Figure 7, in which we plot 
J-'b/T^ vs. l/w. We see that for l/w = 0.4, the flux J^b is suppressed by about a factor of 
2 relative to that for infinite /. We have also verified that the reflection asymmetries at a 
given energy are largely insensitive to the form of the potential beyond the corresponding xq. 
Although our method of including the effects of a finite mean free path is certainly not the 
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final word, it should give a reasonable estimate of the magnitude of the suppression. Note 



in particular that the integrand in ( |3.6| ) is largest for lower energies within the enhanced 
reflection zone. For these energies, Xq is on the symmetric phase side of the wall, and the 
dependence on l/w is not severe. This qualitative explanation is consistent with our result 
that J^B is only suppressed by a factor of 2 for //w = 0.4. 

If / is larger than we have estimated, as for example would be the case for more non- 
relativistic 0's for which D is larger and v,/, is smaller, then the final result could be larger 
than that which we estimate by up to a factor of 2. If, on the other hand, w is larger than 
25/T, the final result will be suppressed both by the effect of wm displayed in Figure 3 
and by the effect of l/w displayed in Figure 7. Of course, w < 25/T would correspondingly 
enhance the final result, for example by a factor of 4 for w = I = 10/T. In the estimates 
that follow, we take J-'is/T^ ~ 5 x 10^^, as appropriate for m = A = fi = T= lOOGeV, 
Am/m = 0.2, w = 25/T, l/w = 0.4, = 0.1 and A(3A(p = 0.5. In the regime in which jFg 
is linear in Vy^, A(3, and Aip, we can write it as 

^r^lxlO-^v^ApA^ . (4.3) 

The dependence on the other parameters is more complicated, as we have discussed and 
illustrated above. 

Next, we give an estimate of the baryon number density carried by particles in the 
region in front of the bubble wall. A (p particle emerging into the symmetric phase from the 
wall begins to diffuse, and the mean distance such particles have travelled from the wall a 



time t after being reflected is x ~ ^2lv^t. In the same time t, the wall itself has moved 
a distance v^t. Defining r as the time that a reflected particle spends in the symmetric 
phase before the wall overtakes it, we find that on average, 

2lvs 

T ^ . 4.4 

Note that this means that on average a reflected particle undergoes 
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scatterings during the time r it spends in the broken phase. Our treatment is only consistent 
for Ai'scatt > 1, and our final result is largest for large Ai'scatt- As mentioned above, although 
non-relativistic 0's have benefits, they also have costs, and we see one here. 

We denote the mean separation between the diffusing particle and the oncoming wall 
during the time the particle is in the broken phase by Ax. This quantity will cancel in the 
final result. Over a range of x in front of the wall given approximately by Ax, there is a net 
baryon number density carried by particles given by 

(4.6) 

We arrive at this estimate by noting that jFg is the baryon number injected into the sym- 
metric phase per unit wall area per unit time and that at any given time, the 0's refiected in 
the previous r are in a region of order Aa; ahead of the wall. Thus, we conclude that every 
point in the universe experiences a baryon density in particles given by ( [4.6| ) for a time 



t* ~ Ax/v,^ (4.7) 

while in the symmetric phase. By this point, it should be becoming clear that ns/s will turn 
out to be largest for small wall velocities. The authors of Refs. [|3>E3] discuss wall velocities 
ranging from 0.02 to 0.4, and other authors have considered velocities as high as 0.9. The 
mechanism we are proposing will be most effective at the lower end of this range, and we 
have been using ~ 0.1 in our canonical example. Note also that we are assuming that 
n'^/T^ is small, and hence we are not including the effect of the baryon number asymmetry 
in the distribution functions used in the calculation of jFg in the previous section. 

The mechanism that generates the baryon number flux carried by (p particles into the 
symmetric phase does not involve any baryon number violation. Therefore, the baryon 
number associated with must be exactly cancelled by a baryon number of opposite sign 
behind the wall. Note that the (p particles, being scalars, are not produced anomalously 
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in electroweak sphaleron processes, and therefore cannot bias such processes. Thus, if this 
were the end of the story, we would have made no progress. However, particles can be 
converted into quarks. This is a model-independent statement, equivalent to the statement 
that the 0's have baryon number 1/3. The rate for 0-quark conversion is, however, model 
dependent. In the supersymmetric case we are using as an example, if the gluino mass Mg is 
less than ^s2i the symmetric phase mass of the heavier squark, then the squarks can decay 
into gluinos and quarks, since the quarks are massless in the symmetric phase. We will not 
assume that the gluinos are this light, however. The first process we consider is scattering off 
a gluino in the thermal bath: squark + gluino quark + gluon. The rate for this process is 
suppressed relative to that for squark-gluon scattering by exp(— M^/T) due to the paucity 
of gluinos in the plasma. It is also suppressed by (T/MgY in the cross-section. The second 
process we consider is gluon + squark — > quark + virtual gluino, where the gluino becomes 
quark + anti-squark or anti-quark + squark. This rate is suppressed relative to ordinary 
quark-gluon scattering by of order as{T/MgY , and by three-body phase space. Defining 1/b 
as the fraction of scatterings incurred by a particle in the symmetric phase that turn the 
(j) into a quark, we estimate that b ~ 200 for Mg ~ 300 GeV, noting again that this estimate 
is quite model dependent. We can now estimate that the baryon number density carried by 
quarks in the region Ax in front of the wall is 

q 6 -''scatt 6 ^ ^6 / a o\ 

b bv^ 

where we have used ( (4.51 ). This estimate is only valid for A^scatt/^ < 1; if A^^scatt/^ > 1 then 
the squark-quark conversion reactions have time to establish chemical equilibrium, and an 
analysis in terms of chemical potentials should be used. For A^scatt/^ ~ !> which is relevant 
for Mg ~ 300 GeV and ~ 0.1, the average is converted to a quark at some point during 
its A^scatt scatterings, which leads us to estimate that the baryon number density in quarks 
is 
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where we have used ( [4 .61) . The density does bias electroweak sphaleron processes. Note 
that a baryon number density in front of the wall cannot be affected by non-perturbative 
QCD processes. This is in contrast to what happens in many other mechanisms. For 
example, if an axial baryon number density (more left handed quarks than right handed 
ones; no net excess of baryons) is generated, this can bias electroweak sphaleron processes 
only if it is not first wiped out by QCD processes. 

The rate per unit volume of baryon number violating sphaleron processes in the sym- 
metric phase is conventionally written as 

r, ~ KatvT^ . (4.10) 

Recent work [^] shows that k is parametrically 0{aw)- The most complete numerical 



simulations done to date |^ suggest k = (29 ± 6)aw In thermal equilibrium, no net 
baryon asymmetry is generated by these processes, since sphaleron and "anti-sphaleron" 
processes occur at the same rate. However, in a region with a nonzero n^, the sphaleron 



processes tend to reduce the baryon number density. In particular |^| 

hsr^-Gnl^^. (4.11) 

Assuming that Ygt* /T^ < 1, which is certainly the case for > 0.005, then the net change 
in ub due to anomalous electroweak processes is 

AnB = -6n%^t* . (4.12) 

Long after the electroweak phase transition, after the universe has re-homogenized, the 
remaining excess baryon number density will be given by An^. Note that in order to end 
up with a baryon density in the universe today, the phase in the mass matrix must be 
such that is negative. A net flux of anti-baryons is injected into the symmetric phase, 
compensated by a baryon flux from the wall into the broken phase. Some fraction of the 
anti-baryon excess in front of the wall is wiped out by sphaleron processes and after the 
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entire universe is swept by the broken phase, a positive baryon number asymmetry persists 
at temperatures below T^. 

Using ( |4.12| ), ( ^4.101 ), (^^91), and (^), and noting that the cosmological entropy 



density at the time of the electroweak phase transition is s = (27r^5f*/45)T^ ~ 55T^ with g* 
the number of degrees of freedom in equihbrium, we estimate that the mechanism we have 
presented yields a baryon to entropy ratioQ 

6 4 ITv^ Tb TON 



Inserting the expression ( ^.3| ) for the baryon number flux jFg, valid in the regime in which 
J^B is linear in Vw, A/5 and A<^, and using Vcf, ~ 0.7 and / ~ 10/T, we obtain 

I^^lxlO-io^A/^Ayp . (4.14) 

s vl 

As we have discussed, this result is obtained for baryon number carrying scalars with sym- 
metric phase masses of 0(Tc) that are non-degenerate by Am/m ~ 0.1 — 0.3. The re- 
sult ( ^.14| ) is sensitive to the wall width w. If we optimistically use ti^ = 10/T instead of 
w = 25/T, the BAU increases by a factor of 4 relative to that of ( ^4.141 ). 



The asymmetry (|4.14|) is at an interesting level. For example, if ~ 0.1, a BAU 



consistent with cosmological observation is obtained for 

ApAip > 0.003, (4.15) 
where we have used k ~ 29aw ~ 1 E^. Taking = 0.1 as we have done is reasonable, but 



wall velocities as low as ~ 0.02 are possible ||4yJT^. Therefore, making the most optimistic 



choice for v^, we find that the scalar baryon number transport mechanism can yield a BAU 
consistent with cosmological observation for Aj3A(p > 1 x 10~^, or even A/SAip > 3 x 10^^ 
if w = 10/T. Making more reasonable choices for w and yields (^4.15|) . We see that even 



^Although our final expression has powers of in the denominator, the derivation relies on ( [4.12 ) 
which is only valid for Tgi* /T^ < 1, and the BAU does not in fact diverge for — > 0. 
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if we use the conservative estimate A/3 ~ 0.01 — 0.03 [^,0, no strong constraints need 
be imposed on Ay9 in order for our implementation of the scalar baryon number transport 
mechanism using top squarks to yield a BAU consistent with cosmological observation. 

We have demonstrated that the scalar baryon number transport mechanism can yield a 
cosmologically interesting BAU, and have done so via a supersymmetric implementation of 
the mechanism. It is therefore interesting to compare our result to those of other authors who 
have studied electroweak baryogenesis in supersymmetric theories. Most recent treatments 
p3|^6| have included one stop with a light broken phase mass, in order to have a strongly 
first order phase transition so that the BAU which is generated is preserved [|^, but they 
have used symmetric phase stop masses which are much larger than Tc. In this regime, 
our mechanism is not effective. These treatments have used particles other than stops and 
have used the charge transport mechanism. Given the new results of Moore [^, they have 
likely underestimated the effects of strong interaction processes which wash out axial baryon 
number. Nevertheless, it seems likely that for masses such that the scalar baryon number 
transport mechanism is efficient, the contribution from charge transport involving particles 
other than stops is comparable to that which we find.0 

The one treatment other than ours that includes reasonably small symmetric phase stop 



masses is the work of Huet and Nelson |T^], and we now attempt a more quantitative 



comparison with their results. For one of their sets of parameters,]] namely m = 150 GeV, 



^AU the mechanisms we discuss in this paper are "non-local" , in the sense that the relevant CP 
violation occurs at the bubble walls, while the relevant B violation occurs away from the bubble 
walls in the symmetric phase. Local mechanisms, in which CP violation acts directly to bias the 
gauge and Higgs dynamics of sphaleron processes, have also been considered ||8-|To|,22| and also 



make a contribution to the BAU. The work of Ref. suggests that for the thin wall case of 
interest in this paper, the local contribution to the BAU is likely small. 

^Note that a completely quantitative comparison between our results and those of Huet and 
Nelson is actually not possible, because they choose to neglect those diagonal terms in of (2.1) 



proportional to and mf. This may be justifiable for m = 150 GeV, and we therefore compare 
our results with those they obtain using this parameter set. For their other set of parameters, 
which has m = 60 GeV, A = // = 50 GeV, their has one negative eigenvalue in the broken 
phase, rendering comparison to results they obtain with these parameters difficult. 
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A = jj, = 50 GeV and Tc = 60 GeV they find a contribution to the BAU due to stops using 
the charge transport mechanism given by ub/s ~ 2.3 x 10^^°f^„(K/K')A/3A(/?. Just as k 
parametrizes the rate for electroweak baryon number violating processes, n' parametrizes the 
rate for strong axial baryon number violating processes. The factor {h/k') is conventionally 



taken to be ~ 1, but the work of Moore |^3|] suggests that it is in fact smaller. For ~ 0.1, 
the BAU ( [4.14| ) generated by the scalar baryon number transport mechanism is a factor 
of 400k' larger than that generated by charge transport involving stops. Although the 



mass matrix used in Ref. [|T^ is not given by (|2.1| ), it seems plausible that when the scalar 



baryon number transport mechanism is efficient, namely for m ~ T, Am/m ~ 0.1 — 0.3, 
it yields the dominant stop contribution to the BAU. Huet and Nelson also consider the 
contribution to the BAU due to charge transport involving particles other than stops and 
find ub/s ~ 6.5 x 10~^v^{k,/ k,')A(3 simpB- This is somewhat smaller than ( [4.14|) , but only by 
a factor of 15k' if ipB ^ Aip. Note, however, that in some models fB is suppressed while 
Aif is not. Nevertheless, a complete treatment of the BAU should include the contribution 
due to charge transport involving particles other than stops. Although perhaps not the 
whole story, scalar baryon number transport yields the dominant contribution to the BAU 
due to stops, and can explain the cosmologically observed value. 

V. OPEN QUESTIONS AND MODEL IMPLEMENTATIONS 

We have given our quantitative conclusions in the final four paragraphs of the previous 
section; the present section is devoted to unresolved questions and to a discussion of possible 
implementations of the scalar baryon number transport mechanism. We have organized our 
presentation of the scalar baryon number transport mechanism in such a way that all the 
parts of the treatment requiring technical improvement were deferred to Section IV, in which 
we restricted ourselves to making estimates. For example, further work is certainly called 
for in the calculation of the diffusion constant D for semi-relativistic strongly interacting 
particles. In addition, our treatment of the suppression due to finite mean free paths can 
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be improved. The reader can surely find other ways to improve the arguments of Section 
IV. Moving beyond the technical, we noted at the end of Section IV that a more complete 
treatment should include the transport of charges other than baryon number carried by 
the scalars of interest. Also, we have neglected thermal contributions to particle masses. 
Since they are of order couplings times Tc, and since the scalars we discuss have masses 
of order in the symmetric phase and somewhat higher in the broken phase, neglecting 
thermal masses seems reasonable in this exploratory treatment of the scalar baryon number 
transport mechanism, and we have left their inclusion to future work. 

The crucial observation that makes the scalar baryon number transport mechanism pos- 
sible is the existence of a broad enhanced reflection zone, a range of incident energies in 
which reflection coefficients and their CP violating asymmetries are large. This arises when 
there are a different number of propagating modes at a given energy on the two sides of the 
wall. A broad enhanced reflection zone may arise in contexts other than that which we have 
considered, for example with scalars that do not carry baryon number. This suggests that 
insights gained from this work may have wider application. 

As noted in the introduction, our main goal in this paper has been to present the scalar 
baryon number transport mechanism, not to address model building issues. We have chosen 
to work within a supersymmetric scenario. Within this context, we now discuss some lessons 
for future model building efforts. It has already been realized that it is desirable for 
one stop to have a zero temperature mass less than the top mass, because this assists in 
making the electroweak phase transition more strongly first order. We now see that it is 
also advantageous to have symmetric phase masses that are of order Tc ~ 100 GeV, and 
that differ by 10-30%. It will be interesting to look for supersymmetric models satisfying 
this criterion. As noted in Section II, non-degeneracy of the appropriate magnitude arises 
in some models due to renormalization group running down from a high energy scale at 



which the stops are degenerate. 

We have called our mechanism scalar baryon number transport rather than stop baryon 
number transport deliberately. The essential features of our mechanism can be implemented 

31 



in other extensions of the standard model involving baryon number carrying scalars, although 
such extensions are perhaps not as well motivated as the supersymmetric scenario. The re- 
cent HERA anomaly |32| may hint at the existence of first generation scalar leptoquarks of 
zero temperature mass ~ 200 GeV. (See, for example, the treatment of Babu et al. [0.) 
A scalar leptoquark is a particle with a Yukawa coupling to a quark and a lepton, which 
therefore carries both lepton and baryon number. Suppose that there are 3 generations of 
leptoquarks diagonally coupled to the 3 generations of quarks and leptons by their Yukawa 
couplings. In a model with two Higgs doublets Hi and couplings of the form HiHj(j)a(pi3, 
where i,j = 1,2 and = 1,2,3, can contribute to the masses of the leptoquarks in the 
broken phase and provide CP violating mixing. If the leptoquark masses receive other con- 
tributions that are nonzero in the symmetric phase, then a mass matrix of the required form 
can arise. (In fact, if contributions beyond tree-level are included, CP violating mixing terms 
can arise even in a theory with a single Higgs field.) The simplest possibility, namely mixing 
between first and second generation leptoquarks, is tightly constrained by bounds arising 
from the non-observation of fiavor changing neutral currents , and probably cannot yield 
large enough off-diagonal terms in the mass matrix to be of interest. However, such bounds 
are absent or much weaker for 02 ~ 03 or 0i — 03 mixing. Implementing the scalar baryon 
number transport mechanism using leptoquarks is therefore possible. It requires symmetric 
phase masses of 0(Tc), but since the zero temperature masses receive additional contribu- 
tions proportional to the Higgs vacuum expectation values, these can still be ~ 200 GeV or 
higher. It is possible to construct leptoquark theories that are consistent with experiment 
in which the scalar baryon number transport mechanism generates a BAU consistent with 
observation. We postpone further investigation, in particular attempts at using the BAU 
to constrain parameters of such theories, until such a time as the experimental evidence 
becomes more compelling. 

Let us hope nature is such that the stop (or leptoquark) spectrum is soon within reach 
of experiment, enabling us to discover whether the mass matrix is such that the observed 
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cosmological baryon asymmetry can be due to scalar baryon number transport. The mech- 
anism is efficient if there are two scalars with symmetric phase masses of order that differ 
by 10-30%, and if the bubble walls during the electroweak phase transition are sufficiently 
slow and sufficiently thin, as we have discussed in our conclusions presented in the previous 
section. Under such circumstances, the baryon number flux produced by reflection of scalars 
with incident energies in the enhanced reflection zone can easily lead to a baryon asymmetry 
of the universe consistent with cosmological observation. 
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APPENDIX A: CALCULATION OF THE REFLECTION COEFFICIENTS 

In this Appendix, we describe a method for numerically evaluating the reflection co- 
efficients used in Section II. We want to find solutions to the time-independent two field 
Klein-Gordon equation (|2.6| ). Solutions to these second order linear ordinary differential 
equations are uniquely determined by specifying four boundary conditions on the fields 
and/or their ffist derivatives. Along with the linearity of the differential equations, this 
implies that the solutions form a linear vector space of complex dimension four. 

As discussed in Section II, because particles incident from the broken phase do not yield 
significant asymmetries, we need only consider the problem of calculating the reflection 
coefficients for 0's incident from the symmetric phase. In order to calculate Ru, for instance, 
we must find a solution that satisfies the following conditions: At large negative x in the 
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symmetric phase, 0i has a right-moving {i.e. incident) component with unit amphtude and 
02 has no right-moving component. There is no restriction on the left-moving plane waves 
in the symmetric phase — their amplitudes determine the reflection coefficients. In the 
enhanced reflection zone, the solutions in the broken phase have one propagating mode, 
which must be purely right-moving, and one non-propagating mode, which must decay 
(rather than grow) exponentially for x ^ oo. 

In the relevant solutions, the propagating modes in the broken phase must be purely right- 
moving {i.e. outgoing), and the non-propagating modes must be exponentially decaying. 
These two broken phase boundary conditions, one for each mode, restrict the space of 
relevant solutions to a two-dimensional subspace of the complete four-dimensional solution 
space. In the basis of interest to us, the two solutions spanning this subspace correspond 
to an incident symmetric phase 0i with no incident 02 and to an incident symmetric phase 
02 with no incident 0i. To find these basis solutions directly requires imposing boundary 
conditions in the symmetric phase. Imposing two boundary conditions at large positive x 
and two at large negative x yields a more time consuming numerical task than imposing 
four boundary conditions at one point. Instead, we proceed as follows. We first find two 
linearly independent solutions satisfying the broken phase boundary conditions, but not the 
symmetric phase boundary conditions. We find each solution by imposing four boundary 
conditions at one point in the broken phase and using the Runge-Kutta algorithm built into 
Mathematica [^]. These two solutions form a basis for the subspace of solutions satisfying 
the broken phase boundary conditions, and we find the basis of interest, namely the solutions 
satisfying the symmetric phase boundary conditions, by taking linear combinations. 

The boundary conditions described above should in general be imposed at spatial infinity. 
We find solutions satisfying boundary conditions at finite x+ in the broken phase and at 
finite a;_ in the symmetric phase. Because we have chosen our our profile function p{x) of 
( p.2[ ) such that the mass matrix does not vary for x < —w/2 and for x > +w/2, we can 
set x+ = w/2 and X- = ~w/2 without loss of accuracy. For a different profile function, for 
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instance one with exponential tails, one would have to choose x_ and x+ far enough out on 
the tails to achieve the desired accuracy. 

In Section IV, we discuss a method for obtaining a crude estimate of the effects of a 
finite mean free path. Instead of imposing boundary conditions at x_ ^ — oo (equivalently 
for our profile function, x_ = — w/2), we impose them at X- = Xq — I. Here, Xq is the point 
where one mode is totally reflected and is found by setting one of the eigenvalues of the 
mass-squared matrix equal to E"^ and I is the mean free path. In order to implement this 
calculation, in the formalism we present below we keep x_ a free parameter. 

We begin by finding two solutions, solutiouo, and solution^j, by matching at x+ to the 
following right-moving solutions: 

solutiouQ,: (f){x) = Ae~'^'"^''ubi 

solution^: = Be~'P'^^''ub2 (Al) 



where pbi = \J E"^ — and pb2 = \J ~ 1A2 ^^e momenta of the normal modes in the 
broken phase at x = x+; ^bi and /i{,2 are the masses of these normal modes, defined as the 
square roots of the eigenvalues of the mass-squared matrix in the broken phase; and Ubi 
and Ub2i the eigenvectors of the broken phase mass-squared matrix, define the broken phase 
normal modes in the symmetric phase cf) = {4>i,4>2) basis. If a mode has real momentum. 



( |Al| ) ensures that it is right moving, since the time dependence of all modes is exp{iEt). If 
a mode is non-propagating, as one mode is in the enhanced reflection zone, its momentum 
should be taken to be negative imaginary to give a decaying exponential and not a growing 
one. Matching to the asymptotic solutions is equivalent to imposing the boundary 

conditions 

solutiouo,: = Ubl] = —ipbiUbi 

solution/3: = Ub2; = -ipb2Ub2 (A2) 



at X = x+. Solutions to the Klein-Gordon equation ( |2.6| ) with these complete one-point 
boundary conditions are easily found. The solutions that satisfy the symmetric phase bound- 
ary conditions are linear combinations of solution^ and solution/3. Note that if w is too large, 
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we will run into difficulty in the enhanced reflection zone. We impose boundary conditions 
at x+, and find solution^ and solution/? by evolving toward smaller x. These solutions in- 
clude one mode that grows exponentially as x is reduced, until x reaches Xq. Therefore, 
if x+ — Xq is too large, the task of finding the solutions that satisfy the symmetric phase 
boundary conditions involves small differences between exponentially large quantities. We 
have found that going beyond w = 40/T requires about 30-digit working precision, and is 
therefore prohibitive. 

In order to obtain the desired linear combinations of solution^ and solution/?, it is nec- 
essary to know the amplitudes of the incident and refiected components of both 0i and 02 
in the symmetric phase for both solution^ and solution/?. At the point x_ in the symmetric 
phase, we denote the eigenvectors of the mass-squared matrix defining the modes 0i and 02 
by Usi and Us2 (orthogonal because M"^ is Hermitian), the masses (square roots of the corre- 



sponding eigenvalues) by fisi and fis2, and the corresponding momenta by Psi = J E"^ — 



and Ps2 = y-^^ ~ /^s2- Then the solutions at x_ will be of the form 

solution^ : 0(x) = A^ie-^^^^'u,! + 5„ie+'P^i"M,i + A„2e-'^^^"M,2 + 5a2e+'^=^"u.2 
solution/? : 0(x) = Apie-'^^'^u.i + Bp^e+'^^'^Usi + Ap2e~'''^'''us2 + ^/jse+^^'^^'u.s (A3) 

where the A's and 5's vary only slowly with x near x_ provided the mass matrix does 
not change much on length scales comparable to the wavelengths of the modes there. This 
condition is identically satisfied for x_ < —w/2, and is very well satisfied at larger values of 
x_ for a wall width 25/T. We now write the amplitudes in ( |A3| ) in terms of and its first 
derivative at x = x_. The amplitude A^j of the incident component of mode j in solution^ 
is given by 

= e+'-^^--^,. • (A4) 
The amplitude Baj of the refiected component of mode j in solutiouo, is given by 



Ba, = e-^^---^:,. ■ ^^^J±^0(x.). (A5) 
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The expressions for solution/3 are analogous. 

The solutions that have incident modes in the symmetric phase that are either purely 
01 or purely 02 can now be constructed from the amplitudes Aaj, Apj, B^j and Bpj. The 
solution with incident 0i and no incident 02 is the linear combination of solutiouo, and 
solution/3 for which the A2 terms cancel, namely 

solutioui: Ap2 solution^ — Aa2 solution^ . (A6) 



Solutioui can be written in the form (|A3|) , with coefficients 



An 


= Ap2Aal 


— Aa2Af3i 




A12 


= Ai32Aa2 


- ^^2^,32 = 




Bu 


= 


— Aa2Bf3l 




B12 


= Af^2Ba2 


— Aa2Bp2 ■ 


(A7) 



These coefficients are the amplitudes of the incident and reflected modes in solutioni. Sim- 
ilarly, the solution with incident 02 and no incident 0i is 



solution2: A^i solution^ — A^i solution/^ (A8) 



with amplitudes 



A21 


= ApiAal 


- A„iv4^i = 




A22 


= Ai3iAa2 


— AalAf32 




B21 


= ApiBal 


— AaiBpi 




B22 


= Aj3iBa2 


— AalBf32- 


(A9) 



We now have all the ingredients necessary to construct the reflection coefficients Rij. (Note 
that as shown in Section II, Rjj = Rji.) 

The reflection coefficient Rij is defined the be the ratio of the reflected 0j current into 
the symmetric phase to the incident 0i current from the symmetric phase. The current 
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represented by a solution (f){x) is i{dx(t>*)(t> — i(t>*{dx4')- Foi" ^ solution Ae*^^* this current 
is 2|ApRe (p). Hence the reflection coefficients are: 








Re psi 


Af32Bal 


— Aa2Bf3i 


2 




All 


Re psi 


A02A,,i 


" A„2A3i 






Bl2 


^ Re ps2 _ 


Ai32Ba2 


— Aa2Bf32 


^ Re Ps2 




All 


Re Psi 


Ag2A,,i 


- A,,2A3l 


Re psi 




B21 


^ Re psi _ 


ApiBal 


— AalBpi 


^ Re Ps\ 




A22 


Re Ps2 


Ai3lAa2 


- A^iA32 


Re Ps2 




B22 


^ Re Ps2 _ 


Ai3lBa2 


— AaiB02 


2 




A22 


Re Ps2 


Af3lAa2 


— AalAfj2 





(AlO) 
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